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Abstract 

, A zero-finding technique for solving nonlinear equations more efficiently than they 

usually are with traditional iterative methods in which the order of convergence is 
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improved is presented. The key idea in deriving this procedure is to compose a 
given iterative method with a modified Newton's method that introduces just one 
evaluation of the function. To carry out this procedure some classical methods with 
different orders of convergence are used to obtain root-finders with higher efficiency 
index. 

Keywords: Newton method; nonlinear equations; iterative methods; order of convergence; effi- 
ciency index. 
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1 Introduction 

No doubt that Newton's method is one of the best root-finding methods for solving 
nonlinear equations. Recent results improving the classical formula at the expense of an 
additional evaluation of the function, an additional evaluation of the first derivative or 
a change in the point of evaluation can be found in the literature on the subject (see 



[HI [TOl |2] and the references therein) . In these works the order of convergence and the 
efficiency index in the neighborhood of a simple root have been improved. 

Using the technique that consists in composing a modification of Newton's method with 
an iterative method we obtain a root-finder for solving nonlinear equations with improved 
order of convergence and efficiency index. The key idea to improve or even double the 
order is to use only one additional evaluation of the function instead of the two evaluations 
needed when applying composition with Newton's method, as it is well-known. 

Currently, IEEE 64-bit floating-point arithmetic is sufficient for the most commonly 
applications in order to obtain the accuracy desired. But, it is increasing the num- 
ber of applications where it is required to use a higher level of numeric precision pQ. 
Namely, evaluating orthogonal polynomials, high-precision solution of ODE's, divergent 



1 



asymptotic series, discrete dynamical systems, experimental mathematics, supernova 
simulations, climate modeling, and nonlinear oscillator theory among others. So, adap- 
tive multi-precision arithmetics facilities are most appropriate in a modern large-scale 
scientific computing environment. 



2 Main result 

Let f(x) = be a nonlinear equation where / : D C R — > R is sufficiently smooth 
in a neighborhood / of a simple root a. Let <f)(x) be an iterative function of order of 
convergence p in I obtained using f(x), f'(x), . . . , f^ p ~ l \x). Hereafter, a technique that 
consists in an iterative method in two steps, is presented. Namely, 

z n = <p(x n ), (1) 

x n +i = z„,- f{z n ) 9g , (2) 

where in (j2J) the factor l/f'(z n ) in the classical Newton's method has been replaced by 
g'q, which is an approximation of the derivative of the inverse function of /. That is, if 
w n = f(z n ), then z n = g(w n ). Therefore, g'{w n ) = l/f'(z n ). As we will see later on, 
this approximation is given by 

/ g{w n )-g{y n ) k -qj k ), . , , k -i 



9 = q 



w n -y n ^ k\ 



+ EV^wk-^ _1 . ( 3 ) 



where y n = f(x n ) and p and q are integers such that p > q > 2. 

Recalling that ([T]) is of pth order of convergence and that we have computed the function 
/ and its derivatives up to order equal to p — 1 at point x n , to analyze the order of the 
two-step iterative method given by (pQ) and ([2]), we state and prove the following main 
result: 

Theorem 1. Let e and E be the errors e n = x n — a and E n = z n — a = K e p + 0(e p+1 ) 
in sequences {x n } and {z n } respectively. Then the order of the iterative method defined 
by d2]) _ (0) is equal to p + q. More precisely, 

\B q+l f'(aY KeP + i\ + 0{eP +q+l ), if p > q, and 



l)«B,+i/'(a) , + A 2 K]Ke 2 P\ + 0(e 2 P +1 ), if p 



f {k \a) g (k) (0) 

where Ak = uj^ry and Bk = wwr k ~ 2 - 

Proof. Putting y instead of y n and w instead of w n , and considering Taylor's develop- 
ments of the functions g(w) and g'(w) in powers of w — y, we obtain 



9 Q {i) (v) ■ <7 (<?+1) (£) 



sh = }_ (" -v)'+ ,_, :„ (« - vr, (i) 

i=0 
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where £ and 77 lie between y and to. From © we get g^ q \y) that after putting it into © 
yields 



^ (i - 1)' 

w-y 

9q + T g> 



(w - yf 



i=0 



(9 + 1)! 



where g' is given in ([3]), as claimed before, and 

9 I ?! (9 + 1)! 



(w - yf 



From © and developing T g in powers of w — y, we have 



g iq+1) (y) 
(9 + 1)! 



+ 0{(w- y y +1 ), 



and from 



9+1 



= a + 5 '(0) y + £ B fe y* + 0{y q+2 ) 



\ fe=2 / 

we have g^ +1 \y) = g'(0) [(q + 1)1 B g+1 + 0(y)}. 
Similarly, writing g(w) = g'(0) [w + B 2 w 2 + <9(u> 3 )] , we get 

g'(w) = g'(0) [l + 2B 2 w + 0{w 2 )] . 
If we develop ([8]) in Taylor's series at point 0, then we obtain 

T q = g'(0)B q+1 (w- y y + O(y0 +1 ) 
= (-l) q g>{V)B q+l y q + 0(y q+l ). 

Substituting © and JTDJ) into © yields 

= g'(0) [l+2B 2 w+ 0(w 2 ) + (-iy +1 B q+1 y* + 0(y q+1 )} . 



Now setting 



y = f(x) = f(a) ( e + Y^A k e k + 0(e^ 2 )\ , 

V A;=2 / 

/(*) = /'(a) (£+ ^ 2J E 2 + 0(i? 3 )), 



we obtain 



5 '(0) [1 + (-iy +1 B q+1 /'(«)« e« + 2 /'(a) £ 2 £ + 0(e^ +1 )] 
g'(0) [1 + /'(«)« e q — 2A 2 E+ 0(e q+1 )] , 



(6) 
(7) 

(8) 



(9) 



(10) 
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where in the last expression we have put f'(a) B 2 = — A 2 . 

Subtracting a from both sides of ([2]) we get e n +i. Assuming that p > q, from the previous 
expression of g' we get 

e n+1 = E — (E + 0(E 2 )) (1 + (-l)**f(a)* B q+l e* + 0(e« +1 )) 
= {-lff'{afB q+1 ^E+ 0{e^E). 
On the other hand, if p = q, then 

e n+1 = E - (E + A 2 E 2 + 0(E 3 )) (l + f'(a)" B q+l e i-2A 2 E + 0{e« +1 )) 

= f'(a) q B q+l e q E+ A 2 E 2 + 0(e q+1 E). 

Replacing E by E = K e p + 0(e p+1 ) the statement follows. □ 

For q = 3 equation ([3|) was used by Kou et al. in [6]. Other contributions related to 
family ([3]) can also be found in [5]. 

Previously, we have set /'(a) -B2 = — A 2 . This relation can be easily proven. From a 
theorem of Jabotinsky [4j, we have 

with the sum is taken over all nonnegative integers such that Y2e=2 ~ Pe = 
and where r = Y2e=2 ^ ne P ro °f °f this theorem can also be found in [9]. The values 
of B q+ \ in terms of Aa 's, for 1 < g < 4, are presented in Table 1. 

Table 1: Values of f'(a) q B q+ i in terms of Aj 
q f'{aYB q+1 

1 -A 2 

2 2^2-^3 

3 -5Al + 5A 2 A 3 -A A 

4 UA 2 l -2lA 2 l A 3 + 6A 2 A A + 3Al-A 5 



We have used the classical definition of efficiency index given in [SJ. That is, EI = 
m l l r ', where m is the local order of convergence of the method and r is the number of 
evaluations of the functions per step. Considering the improvement in the order obtained 
in Theorem 1 the efficiency index is increased considerably. In the case in which the first 
step in the iterative method is of pth. order and there are p evaluations of the functions 
per iteration, the efficiency index is EI = p l l p . By increasing the value of q > 2 and 
applying Theorem 1, we obtain EI = + In Table 2 several values for efficiency 

are given. 



4 



Table 2: Efficiencies 





P 


= 2 


V 


= 3 


V 


= 4 


Method <j>{x) 




« 1.414 




« 1.442 


41/4 , 


« 1.414 


q = 2 


41/3 , 


« 1.587 




« 1.495 


6 l/5 ; 


=s 1.431 


q = 3 








« 1.565 


7l/5 ; 


« 1.476 


q = 4 










8 l/5, 


« 1.516 



3 Some related methods 

In this section, some methods that give the best efficiency indexes for 2 < p < 4 are 
constructed. Notice that in Table 2 the best efficiency index correspond to the case 
when q = p. The expression of the asymptotic constant error for known methods is 
given. Furthermore, we have also computed symbolically and in a different way, the 
asymptotic constant error for the related methods presented and they agree with the 
results obtained using Theorem 1. 



For p = 2, we choose Newton's method as the z = <f>(x) method. If we write 



z n = ip 2 ( x n) = x n - u(x n ) and g' 2 = 2 [y n ,w n ] g - g'{x n ), 

f(x n ) g(w n ) — g(y n ) z n — x n 

where u( Xn ) = ^ and [y n , Wn } 9 = -— = J^J^y then 

Xn+1 — Z n J\Zn) 



or 



f(Zn) ~ f(Xn) f'{x n ) 



,4/ \ _ f(Xn) + f(z n ) f(z n ) 

x n+l — V2\ x n) — z , 



f(Xn) ~ f(z n ) f(x n ) 

Recall that the expression of error in Newton's method is E = A 2 e^ + 03. In the method 
described here the order goes from 2 to 4 and the difference error equation is 

e„+i = A 2 (3A 2 2 -A 3 ) ei + 5 , 
which agrees with the result of Theorem 1 for this particular case. 
• For p = 3, we use Chebyshev's method [2] as the z = 4>(x) method. If we write 



and 



Z n = ^l{x n ) =X n - f 1 + ~ L(x n )^J u( Xn ), 

9s = 3[y n ,w n } g - 2g'(y n ) - - g " \y n ) (w n - y n ) 

= 3 f{Zn] ~ 7(Xn) " 7W) + 27^ U{Zn) ' f[Xn)) 
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f"(x n ) 

where L(x n ) = n u{x n ) , then we have x n+ i = ip 3 (x n ) = z n — f(z n ) g%. Note that 

/ \ x n) 

the error in Chebyshev's method is E = {2A\ - A 3 ) e 3 n + 4 . The error difference 
equation in this improved method is 

e n+ i = (2 A\ - A 3 ) (7 V - 6 A 2 A 3 + A 4 ) e 6 n + O r , 

agreeing again with Theorem 1. 

• For p = 4, the method z = (p(x) considered is Schroeder's method [SJ. Writing 



and 



(•{(.;•„) = ,r„ - ( I + i L(x n ) - i M(a; n ) u(x n ) 2 ) iH.r n ). 

gl = 4 [y„, w n ] 9 - 3g'(y n ) - g"(y n ) (w n - y n ) - g s'"(2/n) - J/n) 
- 1 : ''" ' r " " + ^4 (/(*») - /(*»)) 
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where 

M(x r 



f'(Xn) 



/On) ~ f{x n ) f'{x n ) f'(XnY 

i / n^) 3f{x n f \ , 

6 v. /'(x„) 8 J U[ n) H n)) ' 
/ f"(x n )\ 2 



/"'(x„) . 1 //"(.''„) s 



The error in Schroeder's method is E = (5 A 3 , — 5 A 2 A 3 + ^4) e 4 + O5. The improved 
method presented here is of 8-th order and the error equation is 

e n+1 = (5Al-5A 2 A 3 + A 4 ) (19 A 2 A - 26 A 2 2 A 3 + 7 A 2 A 4 + 3 A 3 2 - A 5 ) e 8 n + 9 , 

agreeing again with it was obtained in Theorem 1. 

Table 3: Test functions, their roots and their initial points 



function 


a 


x 


fi(x) = x 3 - 3x 2 + x — 2 


2.893289 


2.5 


f2(x) = x 3 + cosx — 2 


1.172578 


1.5 


f 3 (x) = 2 sin x + 1 — x 


2.380061 


2.5 


f 4 ( x ) = ( x + l)e- x - 1 


0.557146 


1.0 


/ 5 ( x ) = e ^ 2 +7x-30 _ x 


3.0 


2.94 


f&{x) = e~ x + cos(x) 


1.746140 


1.5 


f 7 {x) = x — 31nx 


1.857184 


2.0 
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Table 4: Iteration number and total number of function evaluations (TNFE) 





Vl 




Vl 


Vl 


vl 


v>t 


VI 


V>2 


VI 


A 0*0 


13 


7 


9 


6 


6 


7 


6 


5 


5 


f2{x) 


13 


7 


8 


6 


5 


7 


•5 


5 


5 


J3{X) 


1 1 
1 1 


D 


Q 
O 





r. 











A 


U(x) 


13 


7 


8 


6 


5 


7 


5 


5 


5 


h(x) 


14 


8 


9 


6 


6 


7 


6 


6 


5 


h(x) 


11 


6 


8 


5 


5 


6 


5 


5 


4 


fr(x) 


12 


6 


8 


6 


5 


6 


5 


5 


4 


Iter 


87 


47 


58 


40 


37 


46 


37 


35 


32 


TNFE 


174 


141 


174 


160 


148 


184 


185 


175 


160 



4 Numerical experiments and comparison 

We have tested the preceding methods with seven functions using the Maple computer 
algebra system. We have computed the root of each function for initial approximation 
xq, and we have defined at each step of the iterative method the length of the floating 
point arithmetic with multi-precision given by 

Digits := p x [-log \e k \ + 2] , 

where p is the order of the method which extends the length of the mantissa of the 
arithmetic, and [x] is the largest integer < x. The iterative method is stopped when 
\ e k\ = \ x k — ct\ < 10~ v , where r] = 3000 and a is the root. If in the last step of any 
iterative method it is necessary to increase the number of digits beyond 3000, then it 
is done. Table 3 shows the expression of the functions tested, the initial approximation 
xq which is the same for all the methods, and the approximation of root a with seven 
significant digits. The functions tested are the same as those presented in [3]. Table 3 
shows the functions; the initial approximation, which is the same for all the methods; 
and the root with seven significant digits. 

In Table 4, for each method and function, the number of iterations needed to compute 
the root to the level of precision described is shown. The notation works as follows: 
Newton's iterative method and the modified method are written as ip2 (p = 2) and ^| 
(p = 2, q = 2). Chebyshev's methods are represented by {p = 3), V| (p = 3, q = 2) 
and ip® (p = 3, q = 3). For Schroeder's method and the modified method we have 
V| (P = 4), V| {P = 4, q = 2), ip\ (p = 4, q = 3) and ipf (p = 4, q = A). In a compact 
way the notation used is ipp +q ■ Notice that the low cost of the iteration functions V4 and 
V|, which show higher efficiency index than the other methods considered. In general the 
results are excellent: the order is maximized and the total number of function evaluations 
is lowest for the iterative methods ift^ and ?/|- 

Finally, we conclude that the methods t/'I and ?/| presented in this paper are competitive 
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with other efficient equation solvers, such as Newton's, Chebyshev's and Schroeder's 
methods (^|, ^3 and ipf respectively). 

5 Concluding remarks 

A technique for accelerating the order of convergence of a given iterative process with 
an additional evaluation of the function is implemented. Furthermore, we have ana- 
lyzed the new schemes obtained from three particular cases: Newton's, Chebyshev's and 
Schroeder's methods. Order of convergence and efficiency index have been improved in 
all these cases. The results have been computationally tested on a set of functions. 

Due to the fact that when the order of convergence of any iterative method is high, we 
need to carry out the computations for testing it with an enlarged mantissa. A multi- 
precision and adaptive floating-point arithmetics with low computing time must be used 
in all the calculations, as we have done in this work. 
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